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We examine the behavior of fermions in the presence of a non-singular thick brane having co- 
dimension 2. It is shown that one can obtain three trapped zero modes which differ from each other 
by having different values of angular momentum with respect to the 2 extra dimensions. These three 
zero modes are located at different points in the extra dimensional space and are interpreted as the 
three generations of fundamental fermions. The angular momentum in the extra dimensions (which 
is not conserved) acts as the family or generation label. This gives a higher dimensional picture for 
the family puzzle. 
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I. INTRODUCTION 

The brane world idea 

0, II H 11 provides an novel 
framework for the unification of gauge fields with general 
relativity and leads to possible new solutions to old prob- 
lems in particle physics which are not resolved within the 
Standard Model: the smallness of the cosmological con- 
stant, the hierarchy problem, the nature of flavor, the 
hierarchy of fermion masses and mixings, etc. 

A key requirement for realizing the brane world idea is 
that the various bulk fields be localized on the brane. The 
brane solutions and different matter localization mecha- 
nisms have been widely investigated in the scientific liter- 
ature [5, 6]. It is especially economical to consider models 
with a purely gravitational localization mechanism, since 
gravity has the unique feature of having a universal cou- 
pling with all matter fields. 

Recently the gravitational trapping of zero modes of 
all bulk fields was demonstrated for the brane solutions 
with an increasing warp factor 0, SJ HI 0- In the 
present paper, using the solution of we show that 
three 4-dimensional fermion generations naturally arise 
on the brane. This gives a purely gravitational mech- 
anism for the origin of three generations of the Stan- 
dard Model fermions from one generation in a higher- 
dimensional theory. The localized fermions are stuck at 
different points in the extra space similar to the model 
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The main point of this paper is to investigate the prop- 
erties of higher dimensional fermions in a 2-dimensional 
curved extra space. In the particular background used 
it is shown that three fermion zero modes occur. These 
three zero modes are taken as a model for three different 
families of fermions. In the present paper we will not 
address the mass hierarchy or the CKM mixing between 
different families. 

We consider a (l+5)-dimensional bulk space-time with 
the brane taken as a 4-dimensional string with a 2- 
dimensional extra space. This solution can be consid- 
ered as a higher dimensional version of the cosmic string 
|12| solution. In 2-dimensional spaces anyonic elementary 
particles, whose angular momentum is not restricted to 
be integer or half integer, are possible 01 ■ Whereas in 
three and higher dimensions all particles must either be 
integer spin bosons or half integer spin fermions, in two 
dimensions particles can have any fractional spin and 
obey fractional statistics. This is because the rotation 
group in two dimensions is U(l), whose representation 
is characterized by an arbitrary real (but not necessarily 
integer) number. A well-known field theoretical example 
of the anyon is a charge- flux composite state in (2+1)- 
electrodynamics |l4|. The anyons imply the existence 
of the new physics in the 2-dimensional space. For ex- 
ample, the anyons are expected to play a crucial role in 
the theory of the fractional quantum Hall effect [15J and 
the behavior of particles in 2-dimensional materials in 
condensed matter Another example is the gravita- 
tional anyon (similar to the case considered here), which 
arises in 2-dimensional spaces with non-trivial topology 
|l7|. This non- trivial topology gives rise to a gravita- 
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tional Aharonov-Bohm effect 01 . A similar effect oc- 
curs for the brane solution considered in this paper and 
it gives the condition determining the angular part of the 
higher dimensional fcrmion wave function. 



II. NON-SINGULAR BRANE SOLUTION 

The action we consider in this article is that of gravity 
in six dimensions 



S = / d 6 Xy/=7j 



Af 4 

—R + A + L 



(1) 



where M , R, A and L are respectively the fundamen- 
tal scale, the scalar curvature, the cosmological constant 
and the Lagrangian of matter fields. All of these physi- 
cal quantities refer to 6-dimensional space-time with the 

signature (H ). 

Variation of the action Q with respect to the 6- 
dimensional metric tensor gAB leads to Einstein's equa- 
tions: 

Rab - \^)abR = j^i (gAB^ + T AB ) , (2) 

where Rab and Tab are the Ricci and the energy- 
momentum tensors respectively. Capital Latin indices 
run over A, B, ... = 0, 1, 2, 3, 5, 6. 

The two extra dimensions are parameterized by a ra- 
dial, r, and an angular, 9 coordinate. We take the brane 
to be located at r = and assume that rotational in- 
variance is preserved. The brane can be visualized as a 
4-dimensional string in six dimensions. In this article for 
the metric in the bulk space-time we consider the follow- 
ing cylindrically symmetric ansatz 0, S @ 



ds 2 = (f> 2 (r) ds 2 (4) - A (r) (dr 2 + r 2 d9 2 ) 



Here 



ds 2 ^ 



9a(3 (%") dx a dx^ 



(3) 



(4) 



is the metric of 4-dimensional space, where the Greek 
indices (a, (3, ... — 0, 1, 2, 3) refer to four dimensions. The 
function A(r) in © must be positive to have space- like 
extra dimensions. 

The 2-dimensional space transverse to the brane usu- 
ally is described by a conical geometry which has a deficit 
angle. Because of this deficit angle one can find nonsin- 
gular solution of Einstein equations only if the energy- 
momentum of the brane is proportional to its induced 
metric pj|. Thus for the stress-energy tensor Tab for 
the brane we take the cylindrically symmetric in the form 



= -9fiuE(r), Tij = -g l3 P(r), T ifJt = , (5) 

where small Latin indices numerate extra coordinates 
— 5,6). In this paper we do not consider the in- 
fluence that matter localized on the brane has on the 



background geometry. Using the ansatz the energy- 
momentum conservation equation gives the relationship 
between source functions 



P' + 4-3- (P - E) = 



(6) 



where the prime = d/dr. 

To solve the equations (J2J we require that the 4- 
dimensional Einstein equations have the ordinary form 
without a cosmological term 



(7) 



Then with the ansdtze © and (J5J the Einstein field equa- 
tions ijSJ) become 
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[P(r)-A], (8) 
[P(r) A] . 



These equations are for the aa, rr, and 99 components 
respectively. 

Subtracting the rr from the 99 equation and multiply- 
ing by <j)/4>' we arrive at 



6' A r 



This equation has the solution 



A(r) 



(9) 



(10) 



where p is an integration constant with units of length. 

System (JSJ, after the insertion of l(TU|) and using JfjJ), 
reduces to only one independent equation. So, taking 
either the rr or 99 component as independent equation, 
we have 



rip <fr 



2rM 4 



71 [P(r)~ A] 



(11) 



In |jj a non-singular solution to with realistic 

source functions was found of the form 



1 + az 2 
l + z 2 



A 



1 



(l + z 2 ) 2 ' 



(12) 



where a > 1 is an integration constant (physically it is 
the value of the warp factor <j> at infinity) and 
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z = - , e 
e 



(13) 
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represent a dimensionless radial coordinate and the width 
of the brane respectively. This solution was constructed 
with the following boundary conditions 



m = a(o) - 1 



o , 



(14) 



where a > 1. These conditions mean that at the ori- 
gin @ has the form of a 6-dimensional Minkowski met- 
ric. Any other values of 4>(0) and A(0) simply correspond 
to a re-scaling of the coordinates. Using the boundary 
conditions (|14|1 at r = the constant p in (|1U|I can be 
expressed by the constants entering the solution (|12|) 



20M 4 



2(a-l) (a-l)A 



(15) 



The solution (|12(l corresponds to the following choice 
of the source functions © 



P(r) = A 
E(r) = A 



4(o + 1) 

5<f> 
3(a+ 1) 



3a 
3n 
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(16) 



At the origin r = 0, where 4>(0) = 1, the source functions 
(|16|l have the value 

P(0) = ^A , £(0) = ^A . (17) 



In the asymptotic region </> (r)\ r:s>e 
functions decrease to 



a > 1 and source 



F(r)l * = ir A<p(0) - 

^WL» £ ^A<£(0). 



(18) 



We want to point out a problem associated with the 
source functions. In general the Einstein equations have 
an infinite number of solutions generated by different 
matter energy-momentum tensors, most of which have 
no clear physical meaning. There is a great freedom in 
the choice of E(r) and P(r); their form is only restricted 
by ©. It is not easy to construct realistic source func- 
tions from fundamental matter fields so that the brane is 
a stable, localized object. The source functions E(r) and 
P(r) given in (|16fl satisfy some physically reasonable con- 
ditions: that are smooth functions which decrease mono- 
tonically as one moves away from r = 0. 



III. DEFICIT ANGLE 



discuss these features as they are important in determin- 
ing the angular part of fermions moving in (|12|) . 

We found that the metric of 6-dimensional space-time 
has the form 



(1 



(1 + z 2 ) 



1 



(dr l + r 2 d6 2 



(19) 



where z — r/e. This solution is similar to the interior 
metric of a cosmic string in harmonic coordinates [20| . 
Performing the coordinate transformation 



z = tan (R/e) 
this metric takes the form 



(20) 



ds 2 = 



R 

cos" I 



R 

II sin" I — 



-dR 2 sin 2 

4 



" S (4) 



(21) 



which is similar to the cosmic string metric found in |l2j | . 

We note that in some 2-dimensional space with cylin- 
drical symmetry (0 < 8 < 2tt) 



dP 



dR 2 + f 2 (R)d9 2 



(22) 



the circumference of a circle of radius R is equal to 
2nf(R), where f(R) is the value of effective radial func- 
tion in the given metric. If the space is flat, we would 
expect that 2irf(R) = 2ttR. On the other hand, in cone- 
like spaces this relation will modified to read 



2tt f(R) = 2ttR 



l-± 
2tt 



(23) 



where 5 is the deficit angle. Solving this equation for 5 
one finds 



<5 = 2tt 



1 - 



W) 

R 



2tt[1 - A(R)} 



where we have introduced the deficit parameter 

f(R) 



A(R) 



R 



(24) 



(25) 



The deficit parameter A for our metric l|21|) has the 
form 



A(R) 



e . 2R 

sm 

2R \ e 



(26) 



From the expression (|20|l we see that effective radial func- 
tion R of our space is finite, since R (r) — > ire/2. 

r — *oo 

From <|25|l we see that close to the origin 



A(r) 



(27) 



The conical geometry of the 2 extra dimensions can 
have interesting global features which are characterized 
by a deficit angle as one traverses closed paths in the ex- 
tra dimensions around the 4D brane. In this section we 



and the deficit angle (|24J) is zero. However, at the core 
radius e the deficit angle is large 



A( e ) 



2 

7T 



0.6 , 



(28) 
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thus the space becomes conical. At infinity 



action can be cast in the form 



A(r) 







(29) 



the deficit angle Ij24(l is 8 = 2tt, this means that space is 
closed. This geometry is closely related to that of the 6D 
metric found by Gibbons and Wiltshire [| , or the conical 
tear drop metric of reference pi) . 

The topology of the 2-dimensional space near the brane 
describes a "distorted cone" or a conical tear drop with 
the brane situated on its tip. At the origin r — ► the 
tip of the cone becomes flat. At infinity r — > oo the 
scale function A — > and determinant of the extra space 
becomes zero. This situation is similar to the case of the 
interior of black hole, in which space is closed upon itself. 



5* = | / d 6 xV=gg AB dA$d B <S> 



K 2 , 
rf<h J 1 



d*aV-5 (4) 0/i* {4) 0" $(4) . ( 32 ) 



where k$ is constant with units of length associated with 
the 4-dimensional scalar field $( 4 ). The integral over the 
extra coordinates in 1|32[) is the same as for the spin-2 case 
above. Thus the integral is finite and the spin-0 fields are 
localized on the brane. 

The 6-dimensional action for the U(l) gauge fields in 
the case of constant bulk components (Aj = const, i = 
5, 6) reduces to the 4-dimensional Maxwell action multi- 
plied an integral over the extra coordinates 



IV. LOCALIZATION OF NON-FERMIONIC 
FIELDS 



Localization of the 4-dimensional spin-2 graviton on 
the brane requires that the integral of the gravitational 
part of the action integral Q over r and 8 for our solu- 
tions l|19fl . to be convergent 



„r, AB r, MN T? J? — 

~99 9 ^AM-TBN — 



(33) 



where ka is some constant with units of length. This 
integral is finite and the gauge fields also are localized on 
the brane. 



M 4 f „ , 

S g = — J dx^^—gR = 



= 2^V( a 3_i) f dJy/ZjQBM , (30) 



V. DIRAC'S EQUATION FOR THE BRANE 

We now turn our attention to the problem of the local- 
ization of the bulk spin- 1/2 fermions on the brane. The 
action integral for the fermion in a curved background 
has the form 



5\j 



% hfr A D B y 



(34) 



where and are respectively the scalar curvature 
and determinant, in four dimensions. The formula for the 
effective Planck scale, which is two times the numerical 
factor in front of the last integral in (|30|l . 



'pi 



-1 



M\aef 



(31) 



is similar to those from the "large" extra dimension mod- 
els • As distinct from these models the effective brane 
width (as) contains the value of warp factor at infinity 
a > 1. This means that for the same fundamental scale 
M the actual brane width e in our case is smaller then 

in a. 

It can be checked that similar to i|30|) integrals of the 
Lagrangian for matter fields also are convergent 0, H 0] • 
This means that zero-modes of matter fields are local- 
ized on the brane with the non-exponential, increasing 
gravitational factor tp. 

Indeed, for the spin-0 fields $, if we assume that they 
are independent of extra coordinates, their 6-dimensional 



where Da denote covariant derivatives, T A are the 6- 
dimensional flat gamma matrices and we have introduced 
the sechsbein hA through the usual definition 



9AB 



h A ^ B 



h B Vab 



(35) 



A, B, ... are local Lorentz indices. 

In six dimensions spinor \1/ has eight components and 
is equivalent to a pair of 4-dimensional Dirac spinors. In 
this paper we use the following representation of the flat 
(8 x 8) gamma-matrices (for the simplicity we shall drop 
tildes on the local Lorentz indexes when no confusion will 



occur 



r r = 



lu 

-7, 

-1 

1 



(36) 



Ta = i 



1 

1 
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where 1 denotes 4-dimensional unit matrix and 7„ , 75 are 
ordinary (4 x 4) gamma-matrices. 

The 6-dimensional massless Dirac equation, which fol- 
lows from the action l|34|). has the form 

(h^T B D^ + h T ~T B D r + h%Y B D e ) *(x A ) = . (37) 

We shall use the following gauge of sechsbein for our back- 
ground metric © 



hP = I 6 B - 6?— S B — 
A r xA' *rVA 



(38) 



Covariant derivatives of spinor field have the forms 

1 



(39) 



do 



From the definition 
2 
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h" M [d M h N N -d N h N M 



ift" 5 (a 



(40) 



the non- vanishing components of the spin-connection for 
the sechsbein (|38|l can be found 



, ,ra _ Xa t 



(41) 



Then Dirac's equation l|37|l takes the form 
> VX V 2r 4A 



1 



-a/A 



V(x A ) = . (42) 



VI. 



ANGULAR DEPENDENCE OF 
WAVE-FUNCTION 



As mentioned in the introduction angular momen- 
tum in two dimensions differs fundamentally from an- 
gular momentum in higher dimensions. In two dimen- 
sions the eigenvalue of the angular momentum opera- 
tor can take values other than integer or half integer. 
In higher dimensions the angular momentum algebra is 
non-commutative, whereas in two dimensions it is a triv- 
ial commutative algebra. We have only one generator 
(S below), which obviously commutes with itself. As a 
result, there is no analogue of the quantization of the 
angular momentum and the only restriction on eigenval- 
ues of the angular operator is the condition to have a 
single-valued wave-function. 



In topologically non-trivial background there generally 
exists the problem of how to find single-valued wave- 
functions. The wave- function ^ of a particle should be 
single valued once we perform a full rotation by the bulk 
angle 9 around the brane. This is non-trivial task since 
geometry described by our solution (|19|l is cone-like with 
a variable deficit angle. 

To find the ^-dependence of ^/(x A ) we need to examine 
the parallel transport of a spinor around the gravitational 
brane solution of l|19|) . The parallel-transported spinor 
at some angle 9 is given in terms of the spinor at 9 = 
by the integration of the condition 



h s °Dg* = h s » ( d e - -luZ ) = , (43) 



where 



and we had introduced a (8 x 8)-matrix which is the gen- 
eralization of the (2 x 2) Pauli matrix 03. 



£ = iT r Tg 



1 

-1 



(45) 



Here 1 denotes 4-dimensional unit matrix. 

Notice that for the ansatz © the 9 component of an- 
gular momentum operator 



u 



do 



(46) 



is not conserved quantity because of factor A(r). If one 
transports a spinor parallel to the circle around the brane 
it will not return to itself but will undergo a rotation 
through the angle, which will be different at the different 
distance r from the brane. So it is not possible in general 
to have a simple separation of the variables r and 9 in 
the wave- function, that will correspond to a single- valued 
spinor on the brane. 

To separate variables r and 9 and find the solution of 
the condition l|43|l we consider small regions of the extra 
space where the deficit parameter (|25f) and spin connec- 
tion (|44() vary slowly and can be considered as constants. 
Under these restrictions the solution of condition l|43l) for 
closed path around the brane is 



^(x A ) = cxp 









iO 


(nA-{ 


-H. 









*(r, x v ) 



(47) 



where n is an integer number. 

Using the intrinsic version of the metric given in i|21|) 
we consider two key regions. The first one is the vicinity 
of the origin r ss (i.e. R as 0), where the metric is flat 



-, A(0) = w(0) = l 



(48) 



Note that in this region the exponent in (|47|l becomes 
16 (n + 5EJ . The term in parentheses is of the form of an 
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orbital angular momentum plus spin angular momentum, 
as is expected near r = where the space is flat. The 
other region we consider is close to the brane surface, for 
example at r — \/3e. At this r the metric (|21[1 describes 
the geometry of exact cone with 



1 

RA 



0.4 



-0.5 



(49) 



In this region we can also separate the variables r and 9 
similar to the case of the cosmic string [l2l 122^ . 

Note that in spite of the fact that the space near the 
origin is flat the spin-connection w is not zero there and a 
spinor parallel transported in a closed circuit around the 
brane picks up an overall minus sign, which corresponds 
to an additional rotation of the spinor about its own axis 
by 2n. This minus sign is a consequence of our choice of 
frame hj? and not a physical effect in case of flat space- 
time. However, for our metric (|19|l at the distance r = e, 
the spin connection (|44|) becomes zero, then changes sign 
and at the distance r = v3e takes the value to = —1/2. 
This changing of spin connection sign is a physical effect 
and corresponds to appearance of anyonic state in our 
model, which can carry fractional angular momentum. 



VII. LOCALIZATION OF FERMIONS 

The condition for the localization of a field on the brane 
is that its "wave function" in the extra dimensions be 
normalizable, or that its action integral over r be conver- 
gent. As remarked before the extra space in our model is 
effectively closed, since R(r — > oo) — > 7re/2, so in general, 
all zero-mode solutions of the bulk fields are normaliz- 
able. However their wave functions may spread rather 
widely in the bulk owing to the lack of an exponential 
warp factor. Thus, in order not to contradict experimen- 
tal constraints such as the charge conservation law, we 
should consider only the modes which are localized inside 
the brane core r < e. This requires the action integral 
over r to be convergent close to the brane surface. 

We approximate the geometry near the origin r ~ 
by a flat metric and near the brane surface r ss e by a 
conical metric. In these regions we can separate the vari- 
ables and the single-valued spinor wave-functions have 
the form (|47|l . Then the Dirac equation Q42|l takes the 
form 



1 



1 



2(j)' 



1 



A' 



-r"d u + -=r r d r + — + — + — - — 



2r 



-a/A 



4A 2r, 
V(r,x v ) = . (50) 



The solution to this equation, if we completely separate 
the 4-dimensional and extra coordinates, is given by 



2 [(u-l)/2+nA]^n^j 
02AV4 ^ z [(v-l)/2-nA]gn( x ») 



(51) 



where ip n (x ,i ) and £"(2;^) are a pair of 4-dimensional 
Dirac spinors and z = r/e is the dimcnsionless coordi- 
nate. The integer number n in (|51|l corresponds to dif- 
ferent possible values of rotational momentum around 
the brane possessing rotational symmetry of the ansatz 
©■ 

We suppose that 4-dimensional spinors ip n (x v ) and 
£ n {x v ) in (|5T|) obey the massless Dirac equations on the 
brane 



^d^ n {x v ) = ^d^ n {x v ) = 



(52) 



Apart from the massless states ij) n (x u ) and ^ n (x L/ ) ) the 
6-dimcnsional Dirac equation can have solutions, which 
correspond to the massive fermions localized on the 
brane. In the single brane models the masses of the 
bounded states are typically of order of the energy scale 
~ 1/e at which the brane (viewed as a topological de- 
fect in higher-dimensional space-time) exists. This means 
that these states are very heavy and we do not consider 
them further here. An interesting supposition is that if 
the energy scale of the extra dimensions is taken to be of 
the order of the electro-weak scale (~ 1 TeV) one might 
envision a scenario where some of the fermions are the 
mass zero modes above and get their mass from some 
Higgs-like mechanism, while other fermions are massive 
modes. This might explain why the mass of the top quark 
is of the electro-weak scale while all other fermions are 
light compared to the electro-weak scale. 

Using the formula for 6-detcrminant for our ansatz @ 



^—g = r^AV^W , (53) 
the action (|34|l of the spin-1/2 fields takes the form 
27rpA 



S* = V I d A x (iB n ^ u d u tp n + 



(54) 



where is some constant with the unit of length and the 
factors B n and C n are proportional to extra-dimensional 
parts of the action integral 

Bn ~ [ dr^z^ +2nA , 



(',, -~ / dr—: 



-2nA 



(55) 



In the expression (|54|l as the angular variable the quan- 
tity 9A was used. These factors represent the behavior 
of the "wave-function" with respect to the extra dimen- 
sions. The condition that we take as indicating a partic- 
ular mode n is localized is for the "wave-function" with 
respect to the extra dimensions to be peaked. In this 
case because of ^-depending exponent in 147|) after in- 
tegration by 9 the mixing between the modes with the 
different angular number n vanishes. Thus in the present 
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approximation, where the A and to are taken as constant 
in different ranges of r, the different generations do not 
mix. A less severe approximation would give the mixings 
between the different generations. 
At the origin r = 



A = oj = <f> = X = l , 
and the factors l|55|) are proportional to 



Bn(0) ' 

C n (0) 



1 



2(1 + n) 
1 



2(l+n) 



y 2(l-n) 



2(1 -n) 
Near the brane surface r — \^3e 



1 



1 



Vx = ^ _ 

<f> 1 + az 2 az 2 ' 
Aw 0.4, -0.5, 

and the factors l|55|) are proportional to 

1 



(58) 



B n (V3e) - 
Cn(VSe) ~ - 



(0.8n- 1.5) 
1 

(0.8n+1.5)* 



(0.8n-1.5) 



■(0.8n+1.5) 



(59) 



To have localization of modes ip n (x"), £ n (x v ) with some 
angular number n inside the brane core r < e we require 
that the factors B n , or C n , should have a maximum near 
the brane i.e. between r w and r w e. This means 
we require -B n (0), C„(0) to be increasing functions (i.e. 
~ z h with b > 0) and B n (v3e), C n (v3e) to be decreasing 



functions (i.e. 



with <i > 0). 



For the case of zero mode n — both factors B n and 
C„ are identical and satisfy the above conditions - near 
r = they are proportional to z 2 and near r = y/3e 
they are proportional to z~ 3 / 2 . Thus two 4-dimensional 



spinors il)°(x v ) and £°(x") (which are undistinguishable 
from 4-dimensional point of view) are localized on the 
brane. For n = 1 only the factor B n satisfies the con- 
ditions for a maximum near the brane (B n (0) ~ z A and 



B n {Vle 



which corresponds to localization of 



^(x"). For the case n = — 1 only the C ra satisfies the con- 
(56) ditions for a maximum near the brane (C„(0) ~ z 4 and 
C„(-\/3e) ~ z -0 - 7 ) which corresponds to localization of 
£ _1 (x l/ ). So totally we have three generation of fermions 
localized inside the brane core. This gives a brane world 
picture for the generation puzzle with the angular mo- 
mentum. 

(57) 



VIII. DISCUSSION 

We have demonstrated that inside the brane core of 
the co-dimension 2 brane solution there naturally arise 
three massless fermionic zero modes. These can be inter- 
preted as a model for the three generations of the Stan- 
dard Model with the angular momentum of the extra 
dimensions acting as the family number. In order to ad- 
dress the fermion mass hierarchy we would need to have 
some mass generating mechanism. One possibility would 
be to introduce a scalar fermion interaction term of the 
form The 4D mass of the fermion is then propor- 

tional to the integral of the extra dimensional part of the 
$(jjr\]/ with respect to the extra coordinates r and 9 as 
in |23 • Another open question of the present approach 
is that the approximation used to separate the 9 and r 
dependence leads to no mixing between the different gen- 
erations i.e. the CKM matrix would be the unit matrix. 
The mixing between the different zero modes and there- 
fore the CKM matrix could be calculated numerically by 
taking into account the variation of u> and A with re- 
spect to r. Both the mass hierarchy and mixing will be 
considered in a future work. 
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